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Abstract 

We study the energy dependence of the saturation momentum in the 
context of the colhnearly improved Leading and Next to Leading BFKL 
I evolution, and in the presence of saturation boundaries. We find that 

the logarithmic derivative of the saturation momentum is varying very 
slowly with Bjorken-j;, and its value is in agreement with the Golec- 
, Biernat and Wiisthoff model in the relevant x region. The scaling form 

^ I of the amplitude for dipole-dipole or dipole-hadron scattering in the 

■ perturbative side of the boundary is given. 

ji^. 1 Introduction 

<D : 

'-^ I The idea of parton saturation was introduced twenty years ago from Gribov, 

Levin and Ryskin as a dual description of unitarity in high energy hard 
^ ■ processes. Since then significant progress has been made in our understand- 

ing of saturation and unitarity 0. Appropriate QCD evolution equations 
that describe high density partonic systems have been derived P, ||, H, 
and saturation is gradually becoming well established theoretically. At the 
same time there has been considerable phenomenological success, as saturation 
based models are used to explain data at RHIC and HERA. The McLerran- 
Venugopalan model J^, ||, § exhibits gluon saturation in a simple and intuitive 
way, and is used in order to describe the early stages of high energy heavy 
ion collisions. In deep inelastic scattering, the model of Golec-Biernat and 
Wiisthoff [0] has been quite successful, as it gives a very good fit to the 
HERA data at small values of x and in a wide region of values. 
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One of the important issues has been the energy (or rapidity) dependence 
of the saturation momentum Qs{Y). How can we determine QsiY) from the 
available QCD evolution equations? To be more specific let's consider the 
scattering of a color dipole of size 1/Q on a dipole or even a hadron of size 
and with relative rapidity Y. Then, considering the forward scattering 
amplitude T{Q, fi,Y), the saturation momentum corresponds to a line in the 

— Y plane, along which T becomes constant and of order that is 

unitarity effects have been reached. 



One should expect BFKL evolution |TT|, (but not its double logarithmic 



limit) to be the relevant dynamics, since this is the one that leads to high 
density partonic systems. For fixed coupling and at leading level this was 
done in [^], and the leading exponential (in Y) behavior of Qg was found. In 
it was reconfirmed and extended to the running coupling case with the 



result log((5^/A^) oc VY. BFKL dynamics, being a linear evolution, suffers 
from the lack of unitarity. Therefore it would seem natural, that the use of a 
more general evolution equation would give the proper description of Qs- This 
was done in [15|, where a numerical study of the Balitsky-Kovchegov equation 
T7| was carried out. 



It would be nice of course to have analytical results from such non-linear 
equations. However, the exact details of the non-linear evolution effects should 
not be important in determining the energy dependence of Qs- As far as the 
system remains on the perturbative side of the saturation line, all the effects 
of the non-linear terms could be described by a simple absorptive boundary. 
Using leading BFKL with either fixed or running coupling, this procedure 



was followed in |]T8[. In the more realistic case of running coupling, negative 
yi/6 corrections to the leading \/Y were found. Even though this correction 
is parametrically small when compared to the leading behavior, it becomes 
important for realistic values of Y, and it tames the growth of Qs which is 
overestimated in the absence of boundaries. 

The goal in this paper is the extension of this procedure at the next to 
leading level. The calculation of the NL correction to the BFKL kernel was 



completed in |T9|, 0. However, this negative correction turned out to be 
larger in magnitude than the leading contribution for reasonable values of 
as = 0.2 -h 0.3. Even worse, when > 0.05, the full kernel has two complex 



saddle points which lead to oscillatory cross sections |2l[. Thus there is no 
hope that one will get sensible and reliable results when using the standard 
NL BFKL dynamics. Perhaps the most appealing cure of these pathologies. 



has been given by Ciafaloni, Colferai and Salam [E2l p3|, E3, E5|. It was rec 



ognized that the large NL corrections emerge from the collinearly enhanced 
physical contributions. A method, the cu-expansion, was then developed to 
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resum collinear effects at all orders, and the resulting improved BFKL equa- 
tion was consistent with renormalization group constraints and the DGLAP 
equations |2^, ^ by construction. The kernel is positive in a much larger 
region which includes the experimentally accessible one. One of the outcomes 
of this method was a very reasonable value for the hard Pomeron intercept 
uop. Therefore we expect the application of the tu-expansion in the saturation 
problem, to give reliable results for the energy dependence of Qs and for the 
form of the scattering amplitude in the perturbative side of the saturation line. 

In Section 2 we give a brief review of the cu-expansion approach, and the 
full form of the NL in uj kernel that we use is given in the Appendix. In 
Section 3 we define a critical line QciX) (not the saturation line) which belongs 
in the perturbative region and which is a line of almost, but not exactly, 
constant amplitude; there is no exponential (in Y) dependence. In Section 
4, by considering perturbations around this line, we approximately solve the 
RG improved BFKL equation under certain but valid assumptions. In Section 
5 we introduce the boundary, which essentially cuts out all paths that enter 
the saturation region Q ^ Qs during the evolution. We remove almost all 
the freedom left in the solution given in Section 4, by fixing the amplitude 
to be constant on the saturation line. The F-dependence of Qs is given and 
a multiplicative constant is free. However, the logarithmic derivative = 
dlog{Q'j./A'^)/dY is well defined. In Section 6 we give our results for As and 
comment on the value of Qs- For dipole-hadron scattering and at NL level in 
u we find Xs{Y = 5 9) = 0.30 -r- 0.29 and that A^ is decreasing very slowly 
with Y. We also comment on the parametric form of the boundary correction 
to As, which is oc [oisiQl)]^^^ and well defined for any value of Y, since Ql is 
increasing during the evolution. For a suitable but reasonable choice of our 
free multiplicative constant in Qs, we can fix its value to be, for example, 
Qs{x ~ 3 X 10~^) ~ 1 GeV. In Section 7 we give the scaling form of the 
amplitude on the perturbative side of the saturation line in terms of Q^ IQ\- 
This is the same (but with a different Qs) as the one obtained in leading 



BFKL with either fixed or running coupling [|TB|, and therefore we discuss its 
generality. In Section 8 we check our approximations and estimate the errors 
which are small and under control, and finally in Section 9 we conclude. 



2 Review of the RG improved BFKL kernel 

In this section we describe as briefiy as possible the collinear improvement of 
the BFKL equation. The full derivation and the details can be found in the 



papers of Ciafaloni, Colferai and Salam p2|, |23|, |23, |23] . The BFKL equation 
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may be written as 



dT 

— = a,K®T, (1) 

with as = O-sNc/ii., and where we can imagine, for example, that T corresponds 
to the forward scattering amphtude for a dipole of size l/Q on a dipole of size 
1/ H- The eigenfunctions of the leading kernel are (Q^/ with correspond- 



ing eigenvalues 12 



Xo(7) = 2^(l)-^(7)-^(l-7), (2) 
where ?/'(7) = d\ogV{'~^) / d'~^. The pole structure of Xo(7) is 

Xo(7) = - + 7^ + finite. (3) 

7 1-7 

The NL correction to the total characteristic function x{l) is o.sXi{l)i where 
Xi{l) is a known function [jl9], |20| and is analytic in the region < 7 < 1. For 
simplicity let us consider its pole structure when Nf = (we will keep Nj 7^ 
later), which is 

where 

127r ' ^ ' 

and Agg{uj) is the Mellin transform of the non-pole part of the gluon-gluon 
splitting function (the pole 1/z has been taken into account in the leading 
kernel), that is 

Agg{uj)= [ dzz^Pgg{z) --. (6) 







When Nf = one has Agg{0) = —b = —11/12 < 0, and one can see in (|) 
the large negative corrections arising from both the double and triple poles. 
It was observed that the poles at 7 = (7 = 1) arise when one considers the 
collinear (anti-collinear) limit ^ iJ? {Q^ -C /x^). 

Since a collinear small-x branching resulted in a pole I/7 in the leading 
kernel, a sequence of a small-a; and a non-small-x branching, as can be seen 
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from (P) will result in a pole 



^-^..(-) - ^«sXo(7)A..(0) - (7) 

where we have truncated the series in at NL level by using uj = agXil:^)- 
This is the double pole at 7 = in Eq.(^). Similarly we obtain the double 
pole at 7 = 1, where there is an extra contribution coming from the running 
of the coupling, if the extracted coupling in Eq.(|l|) is evaluated at Q^. 
However, if we simply do not truncate the result in (^) we should obtain the all 
orders coUinear contribution to the BFKL kernel. Thus, including the leading 
kernel poles, the total eigenvalue will become 

X(7, ^) = + 7^^^^ + finite, (8) 

7 1-7 

which is cu-dependent. This form will be further modified after the following 
discussion, since it remains to interpret the triple poles in Eq.(^). 
Let us write the amplitude T as 



^ ~ ' 27rt 2m ""^^ 



ujY - {1- 7) log % 



(9) 



where Y = log(s/so) and with s the total energy of the scattering, sq is an 
energy scale which depends on the scales Q and ; sq = so(Q, /i). At leading 
level its choice is irrelevant. At NL level one chooses sq to be the natural scale 
of the problem under consideration. In the calculation of cjp a suitable choice 
is So = Qn- In the coUinear (anti-collinear) limit the right choice is sq = 
{sq = n"^). In this limit Y = log(l/a;) = log(s/Q^), where x is now the familiar 
Bjorken variable. As a consequence, the eigenvalue of the NL kernel depends 
on So, and the one given in Eq.(^) is for the symmetric choice So = QyU. From 
(I) we can see that we can switch to Sq = (sq = yU^) by a simple shift 
7 ^ 7 - f (7 ^ 7 + f ), since 

^log^ - (1 -7) log ^ ^ ^log^ - (1 -7) log (10) 
The change in the eigenvalue can be obtained from the following iteration 



^ = asX(7 - ^) ^ (^sXil) - «s^x'(7) (^sXil) - ^tt^Xo(7)Xo(7) 



^^(7) = ^ - 2(l-\)3 + fi^^^^' ^1^) 
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where the result is again truncated at NL level. When we add this correction 
to (^) the triple pole at 7 = cancels, as it should in the coUinear limit, in 
order to have consistency with the RG equations. This is the origin of the 
triple poles in the NL BFKL kernel. Now one can see that Eq.(^ has to be 
substituted by 

when So = Q/i, so that 

7 1 — •y + u 

when So = Q^. 

It remains to combine the behavior (p!2D or (|1^) with the known L and 
NL BFKL kernels. As explained in the next section, the correct choice in 
the saturation problem is sq = Q^. In the cj-expansion approach, where one 
does not truncate at NL in as in Eqs.(0) and (|TT1) , and u appears as the 
expansion variable, the leading kernel is 

Xo(7, ^) = 2^(1) - V^(7) - - 7 + (14) 

or 

Xoil, ^) = Xoil) - -r— + — ^ — , (15) 
1 — 7 1 — 7 + u; 

where Xo(7) is the one given in Eq.(|^). Both descriptions are equally good. 
Then we consider the NL in uj kernel. We need to add to add in Xi{l) ^ term 
-Xo(7)Xo(7)/2 since Sq = (see Eq.(|n|)). After subtracting the NL BFKL 
corrections that are already included in (0) or (p!5D to avoid double counting, 
we replace all remaining poles at 7 = 1 with poles at 7 = 1 + (-01, making at the 
same time the substitution Ac,c,(0) Agg(co') and we obtain xi(7,co'). Then 
the NL in lo eigenvalue is 

x{l,^) = Xo{l,^)+^—i (16) 



which has the proper limit given in Eq.(|TB|) and includes the NL BFKL cor- 
rections. The details in constructing xi(7,ti;) are given in the appendix. 

Our freedom in using either ([T^) or ([TsD for Xq{1i^) (and similarly for 



Xi(7,u;)) emerges from our freedom in how to shift the poles; one can always 
add a finite part. However, it was shown that any resulting difference is of 
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order 0{uj'^). This coUinear resummation independence was checked in the 
value oi Up where the results differ not more than 2 — 3%. We expect this to 
be true in the saturation problem too, since the two kernels are very close to 
each other for all allowed values of 7 and u. For practical reasons we shall use 
a simple pole shift as in (0), since in this case the cj-dependence becomes as 
simple as possible. Notice that because of this a;-dependence, the operator K 
in Eq.(0) contains 8/ BY derivatives. 



3 Evolution along the critical line 

As mentioned earlier we consider the forward amplitude T{Q, /i, Y) for dipole- 
dipole or dipole-hadron scattering at NL level in uj. Since the kernel eigenvalues 
are x{li^)i the evolution equation may be written in an operator form 

X U + T^>^ (17) 



dY hp'^y dp' dY ^ 

Here we have defined the logarithmic variable p = log(Q^/A^), with A the usual 
QCD parameter. The rapidity Y is given by F = log(l/x), with x = Q^/s 
and s the total energy of the scattering. It is obvious that we have chosen an 
asymmetric energy scale Sq = Q"^. The reason is that we are going to evolve 
the system, following a particular line in the — Y plane along which is 
increasing (we should emphasize though that the system is not in the coUinear 
regime). Furthermore, this is the standard expression of Bjorken-a; which is 
also used in the HERA data. The factor 1/bp is the running coupling a^, with 
b from (^. 

We change variable from p to rj according to 



rj = p-p,{Y), (18) 



so that 



d d . d d d 

The dot represents a derivative with respect to Y . The function PciY) will be 
defined soon to be a hne of almost, but not quite, constant amplitude; there 
will be no exponential dependence in y. 77 corresponds to fluctuations around 
this line. If we also expand the running coupling as 



hp bpc hpi \p 



(20) 
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Eg ■ (|I7D becomes 



It is convenient to view the eigenvalue xily ^) ^ function of 7 and A, where 
A is given by 

u; = A(l-7). (22) 

Then, if 

0(7,A)=x(7,A(l-7)), (23) 

Eq.(pTD may be written as 




(T/a,). (24) 



We drop the d/dY in the second argument of the function. This is a vahd 
approximation that will be checked at the end. Thus the second argument 
will be pc- Even though pc is expected to be small, we do not expand around 
0, since such an expansion would lead us to the usual BFKL evolution, where 
the coefficient in the expansion would be large at NL level. We do expand 
only with respect to the first argument, in a "diffusion approximation" much 



like that introduced in |29], around the point 7c (^) to obtain 



d 



1-^ 



0(7c, pc) 



drj 



+ 0l'"(7c, Pc){l + T-~lc]+ T^r'^lc, Pc) 1 + - 7, 



d 



dr] 



(T/a,), (25) 



where 0"*'" is the (m, n) derivative of (/)(7, A). 
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At this stage we need to determine the F-dependence of 7c and pc- We 
choose them, so that T is as constant as possible. This will happen if the 
coefficients of the constant and d/dr] terms in (pSD vanish. Thus 

0(7c, Pc) + (1 - 7c)0''°(7c, Pc) - Whpc = 0, (26) 

and 

Wc + </'''°(7c,Pc) = 0. (27) 

Because of the last term in (p6| ) the coefficient of the constant term in ( pSf ) 
does not vanish. Nevertheless, its appearance is required in order to cancel 
F-prefactors that will appear in the amplitude, by a suitable choice of the pure 
number j3. At the end we will find (3 = 1. Eqs.(p6D and (^) cannot be solved 
in general, except in the simplified cases of collinear models where one 
keeps only the pole parts of the characteristic function x{li^)- These models 
could be useful in studying the form of QsiX) (which is not so simple even in 
this case), but its value is not very close with the one obtained when we use 
the full kernel. However, one can expand around the point 70, which is the 
solution to 



Xo(7o) + (l-7o)Xo(7o) = 0, (28) 
and with numerical value 70 = 0.372... Notice that (03) is the asymptotic form 



of (^6]) when Y is very large, since from (^) we expect the asymptotic behavior 



Pc oc 1/yY. Therefore, a first order expansion around the asymptotic value 
7o, simplifies our equations to 

, -0(70, Pc) - (1 - 7o)0^'°(7o, Pc) + 2(3hpc , . 

(1 -7o)0^'"(7o,Pc) 

and 

PcPc = #^ - (30) 
6(1 -70) 1 - 7o 

The numerical value of 0^'''(7o,Pc) is not too far from x'oilo) = 48.5... The 
appearance of such a big number in the denominator of (^) suggests that the 
expansion around 70 should be valid even in the intermediate region Y = 6-r8, 
which is experimentally accessible. 

One can see that for very large Y, the scattering can be described just 
by leading BFKL evolution with running coupling, since pc goes to 0. This 
is easily understood. During the evolution pc is increasing and therefore the 
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value of the running coupling is decreasing. But this approach to leading 
BFKL dynamics will not be realized practically, until Y reaches a very large 
and inaccessible value. 

For what follows in the next two sections, we will not need the explicit 
form of Pc(y) and ■jciX)- We return to Eqs. (]29|) and ( pOD in Section 6. 



4 Solving the BFKL equation 

In this section we give an approximate solution of the RG improved BFKL 
equation in the vicinity of the line pdY)- Making use of (l26| ) and (l27|), the 
amplitude equation (^) becomes 



d 



oY V PcJ PcOr] 



1] d 



-^0n7c,Pc)(l-^ 

Pc 



1 + — 

dr] 



{T/as 



0. 



(31) 



When rj <^ Pc one can drop the corresponding term in the first parenthesis of 
the last term. Now let's isolate the leading exponential, in the scaling variable 
7], behavior of the amplitude and write T as 

T = «,exp[-(l - 7,)r7 + 2/51ogp,]7/^(r/, Y). (32) 

We have neglected any target dependent (/i-dependent) prefactors in T, which 
are proportional to 1/p^. We easily find that ip satisfies 



-0^'°(7e,Pe)^ + 0^'°(7c,Pe)^^ 

opc pc or] 

-^0''°(7c,Pc)|^ + micPc) + Vc7j^}^ = 0, (33) 

which becomes now our basic equation. Since there is no explicit Y dependence 
in (|33|) , we have decided to view the function ^ as a function of pc instead of 
Y, ip = ip{ri,pc.), by means of the transformation 

d d 

and we have simplified as much as possible the coefficients by using Eqs.(p^) 
and (13). 
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Formally these coefficients should be expressed in terms of pc, which can be 
done, in principle, when one uses Eqs. (PI?D and (PDj). Practically it is easier to 
have them as functions of Ac = Pc- In any case, even though these are compli- 
cated functions, we will need only their asymptotic values in the approximate 
solution of (|33|). 

It is convenient to make a final change of variables from (?7,Pc) to 
where 



and 



Dp. 



1/3 ■ 



(35) 



t = 3Do{l--fo)p 



1/3 
c ' 



with 



D 



<P''\lc,Pc) 



1/3 



and 



1/3 



(36) 



(37) 



(38) 



.2Xo(7o). 

Clearly Dq = 1.99... is the asymptotic value of D. After some straightforward 
algebra Eq.(P3|) becomes 



tlj = 0. 



(39) 



While the coefficients of the last two terms have been totally simplified, the 
ones hi{t) and h2{t) of the first two terms are extremely complicated functions, 
given by 



hi{t) 



and 



h2{t) = hit) 



Do(l-7o)0i'Q(7e, Ac) 
/^0(7c,Ac) 

3Z)'(Ae)pc(Ac)" 



1 + 



(40) 



(41) 



4I^P^(Ae) 

As said before, hi and /i2 have been expressed in terms of Ac = Pc- The prime 
stands for a derivative with respect to Ac, and 7c and pc are given as functions 
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of Ac from Eqs.(p9D and (|30|) . 7c which is contained in the definition of can 
also be easily expressed as a function of Ac through a chain rule differentiation 
giving 7c = 7c(Ac)Ac/ Pc(Ac)- In the asymptotic region of large t, or equivalently 
of small Ac, hi and /i2 approach 1, 

lim hi{t) = lim h2{t) = 1, (42) 

with corrections of order 0{l/t^), which turn out to be small not only para- 
metrically, but even numerically for intermediate values of Y. 

It appears difficult to find an exact solution of Eq. (^) even if we set hi = 
/i2 = 1- However, we can find an approximate solution of the form 



= Ai(^-a)exp 



-jh2{t) 



a 



dt 
'hit) 



dt h2{t) 
thiit) 



where a is a constant to be determined and Ai is the Airy function, 
expression satisfies Eq . (|39|) up to a term R{^,t)ip, where 



e 
t^ 



Ahl + hih2 - thi- 



dh2 
dt 



5e 

t^ 



O 



(43) 



This 



(44) 



It is important that this remainder contains no terms of the order f{^)/t, 
as those terms would be important when one integrates (|39D over t. The 
approximate solution given, should be the one that matches the exact solution 
in the intermediate and asymptotic region. It contains the Airy function pT 



and the diffusion type factor exp(— ^^/t), which appear in general when one 
considers BFKL evolution with running coupling. We also mention that at 
this point we loose track of the overall constant which appears in front of the 
amplitude, as a consequence of our inability to solve (El) exactly. 



5 The saturation boundary 



Up to target dependent factors, the amplitude T is given by Eq . (p2D with ip 
as found in Eq. (^31) . It contains two parameters j3 and cr, which are uniquely 
determined from the saturation boundary conditions in the asymptotic region. 
To that end we express T in terms of the variables (^, t) in the asymptotic 
region, that is 



T = Ai(,^ — cr) exp 



3 



at+ (6/3-5)logt 



(45) 
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When ^ becomes negative, the amphtude grows strongly, it reaches a maximum 
value and then it decreases rapidly to vanish at ^ = ^i + o", where = —2.33... 
is the first zero of the Airy function. A decreasing amplitude is of course an 
unphysical situation, which simply means that we loose our ability to describe 
the evolution of the system by using hnear equations. Thus we need to put a 
saturation boundary before T hits its maximum value, along which we fix T 
to be constant, in order to illustrate the effects of the non-linear terms in the 
simplest possible way. Now writing ^ as 



^ = ^i + a + 5^, 



(46) 



we have 



cr 



Ai'(6 



(47) 



and we keep only the linear term as we will shortly see that 6C, is small in the 
region of interest. Then Eq. (^5[) becomes (dropping again unimportant for our 
purposes multiplicative constants) 



exp 



<«i±i^ + (6^-5)logi-^p + logi5 



(4J 



The maximum of this expression occurs when 6C, = 3/t, thus we assume a 
saturation boundary at 



3(1 + c) 
t 



(49) 



where the amplitude is a finite fraction of its maximum value, c is a positive 
constant of order 0{1), that we are not able to determine. On the boundary 
the amplitude is 



exp 



(6 + 4a)t 



+ 6(/5- l)logt 



Clearly, we need to choose our two parameters to be 



6 

' 4 



(50) 



(51) 



so that T is well defined and t-independent on the boundary. 
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Now we are finally in a position to give an expression for the saturation 
momentum. Using the results stated in Eqs.(^DD and (^) and recalling the 
definitions ([T8|) , (|35D , (|36D and (^) we obtain 



P. = Pc + ^/^oPr + ^ 



which gives 



Qs = A exp 



7o 



1 + c 



(52) 



(53) 



2(1 -7o) J ■ 

where we have replaced D by its asymptotic value Dq. It is obvious that in our 
approach there is a free multiplicative constant. However, for the logarithmic 
derivative of the saturation momentum, Xg = Ps, there is no such ambiguity. In 
Section 8 we will discuss possible uncertainties in Eqs. (|5^) and (|53D, as far as 
the energy dependence is considered, which turn out to be parametrically and 
numerically small. In order to make use of (|52D or (|53|) we have to determine 
Pc{y), a task that we now turn to. 



6 Results for 

Given a kernel, we need to solve Eq.(^) to find the F-dependence of pc- In 
general it appears easier to solve (when a solution is possible) for Ac = Pc- 
Considering pc as a function of Ac and differentiating with respect to Y using 
the chain rule, we have 

A. = p;(Ac)^, (54) 
and by separating variables we obtain 

dX,P^ = Y, (55) 



Ac(0) A, 



where 



Pc(Ac) = #^ - (56) 
6(l-7o)Ac l-7o 

as determined by Eq. (^) . 

It is clear that one needs an initial condition specified in order to proceed. In 
leading BFKL with fixed coupling, the initial condition is totally irrelevant 
for Ac since the leading behavior of pc{Y) in that case is purely linear. In 
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the running coupling case, either in BFKL or in cj-expansion, there is a an 
integration constant Yq or Ac(0) as it appears in (p3D, which affects the value of 
Xc{Y). However, as we discuss shortly after and in Section 8, small changes in 
the initial conditions will not change our results significantly and furthermore 
this change is predictable. The most natural initial condition to adopt is 

p,(r = 0) = log^, (57) 

with of the order of the target size. Then it is straightforward to find 
Ac(0) from (^). Now after determining Ac and/or p^, one can finally get the 
logarithmic derivative of the saturation momentum from the differentiation of 
(^), which gives 

K = Xc+^Xc. (58) 

This is our main result for A^. In the following we examine the cases of Leading 
BFKL, Leading in u) and NL in uj dynamics. 

6.1 Leading BFKL - review 

In this case it is straightforward to solve directly for Pc(^) using Eq. (|30|) . We 
find 



fv\ / 2xo(7o)^ , 



pM 



1 - 7o 



1 - 7o' 



(59) 



from which it is trivial to obtain Xc{Y). Then (^) provides us with an ana- 
lytical expression for \s{Y). 

It is very useful to examine two particular limits of Eq. (|59|) , in order to get a 
better understanding of the initial conditions. 

(i) Let's imagine first that the target is a very small dipole, such that the 
second term of the square root in (|59|) dominates, that is 

m » ^^Y. (60) 
Then one can easily see that (^) reduces to the exact fixed coupling result 

Pc(r) = pc(o) + (61) 

b{l - 7o)Pc(0) 
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with Pc(0) as given in (|57|) and the natural identification l/bpc{0) — > cxg- This 
is expected since ( pOD simply states the fact that, the rapidity is not large 
enough (as 1/al parametrically) and therefore the NL and the running cou- 
pling corrections are unimportant. 

(ii) Now let's consider that Y is large enough so that the first term of the 
square root dominates, that is 

Then (BPf) reduces to 



In this case, where the running coupling dynamics is necessary, it is clear that 
the initial conditions are lost in the region that ( p^j ) is satisfied. Parametrically, 
the uncertainty in Xc{Y), due to the initial conditions, vanishes as l/F^/^. 

Thus it becomes crucial to specify the initial condition, e.g. whether the 
target is a small dipole or a hadron, and the region of Y values that one is 
interested in, before making any kind of expansion in (|59D . Even though we 
would like to imagine in general that Y is large, let's consider a realistic case 
where the target is of size /i ~ IGeV and the rapidity is in the region Y = 6-r8. 
It is easy to check that the two terms in (^) are of comparable size. Therefore 
we use this equation without making any approximations. But there are two 
important points emerging from the above analysis, which will also be true 
in the case of the u expansion. The first point is that a running coupling 
description has to be considered in a correct (at least conceptually) treatment 
of the saturation problem. The second point is that the initial condition is not 
lost in our case, that is the result depends on Pc(0), but any small changes in 
the initial conditions are almost unobservable since they parametrically vanish 
as l/pl^"^ when one considers Ac. 

6.2 Leading in u 

We consider the kernel given in Eq. ([T5|) . The function is 
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The pc equation, Eg . (pOf ) , becomes 



. _ Xo(7o) pc 2pc , . 

Kl-7o)^(l + Pc) 1-70- ^ ^ 

Even though pc is expected to be small, we prefer not to do any kind of 
expansion in the second term of the right-hand side, in order to be consistent 
with the cu-expansion approach. If we expand this term, we loose at the same 
time the correct coUinear behavior as explained in Section 2. We also notice 
the negative sign of this term, which is expected, since the leading in u kernel 
contains a part of the NL BFKL kernel. Now, following the general strategy 
presented in the beginning of this section, we can do the integral in (^) to 
obtain 



""'^^'^ ' ^ogi±^--— )-[A.-A.(0)]=y. (66) 



26(l-7o)A2 6(1 -7o)' V A, 1 + A 

Even though we cannot analytically solve this equation for Ac (we would like 
to avoid an iterative procedure for reasons that we explained before), this is a 
transcendental equation that can be solved numerically for a given Y. Then 
we proceed to determine pc from (^) and A^ from (^Sp. 

We mention here that if we use the "equivalent" kernel given in (0), we 
will not be able to reach a result as simple as the one stated in (pBD, as it 
appears difficult to integrate Eq. (|55D . However the two kernels should give the 
same results. 



6.3 Next to leading in cj 

The kernel in this case is the one in (0), where xiil^^) is a known function 



that is constructed in the Appendix following the procedure of Refs.p2l p3l 24, 



25| . Here it is not possible to give a relatively simple expression as Eq.(|66D in 
the previous subsection. The reason is that xi(7,ti;) contains a piece Agg{u), 
the Mellin transform of the gluon splitting function. This has a term ijj{2 + uj) 
and Eq. (|55D cannot be integrated . In the region of interest the expected value 
of a; is c<j = Ac(l — 7c), which can be found (a posteriori) to be around 0.25. 
Thus we can expand Agg^u), or just its 4'{2 + u) part, around u = 0. Then, 
depending on the order of this Taylor expansion, we find a more and more 
complicated transcendental equation (even at 0th order) containing fractions 
of polynomials and logarithms of Ac. This equation is the analogous of (|66|) , 
that we found in the leading in uj case. There is no reason to present such 
a complicated expression here, but it is again an algebraic equation which 
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can be solved numerically. We find, for a given Y, that Ac reaches rapidly a 
fixed value as we increase the order of the expansion mentioned above; even 
if we keep just the 0th order term, the result is very close to the exact value, 
something which is expected since we have approximated only a small piece of 
the NL in u kernel. In our results for As, that appear in the table of the next 
subsection, we have kept enough terms in this expansion so that the obtained 
value of Ac is practically the exact one. A^ is obtained in the way we have 
already explained in the previous subsections. 

There is another, but clearly approximate, iterative procedure, that we can 
follow in order to obtain a simple equation, which nevertheless gives the same 
results for the NL in uj corrections. We write the full kernel as 

Xil,^) = Xoil,^^) (67) 

and treat the term a; 5 as a perturbation in the following sense. For a given 
Y we solve the leading in u problem. Then for the resulting value of Ac 
we determine the number 6 appearing in (|67|) from 6 = xiilc ^^c) / Xoilc ^c) , 
where Uc = Ac(l — 7c) as dictated from (j2^). Then we insert this value of 6 in 
(|67|) and solve the new problem, which is of equal difficulty as the leading in u 
problem. We obtain a new value of Ac and we can repeat the procedure. This 
converges rapidly and after a couple of iterations 6 and Ac become constant. 
Then we can repeat for different Y. We have decided to extract explicitly a 
factor of UJ in (|67D , instead of treating all the NL correction as a constant, 
since this factor characterizes the very leading behavior of the NL term. 

6.4 Summary of results and discussion 

We assume an initial condition Pc(0) = 2 log 5, so that the target is of size 
yU~^ ^ 1 GeV~^, when A ^ 200 MeV. We set the number of active flavors to be 
Nf = 3. In the table we summarize the values of As, obtained for the different 
kernels considered and for various values of Y. 
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Values of A, = d\og{Ql/ K^) / dY 


Y 


BFKL 


L in cij 


NL in uj 




5 


0.433 


0.366 


0.300 


-0.20 


6 


0.419 


0.358 


0.297 


-0.20 


7 


0.406 


0.350 


0.294 


-0.19 


8 


0.393 


0.343 


0.291 


-0.19 


9 


0.382 


0.336 


0.288 


-0.18 












25 


0.278 


0.258 


0.238 


-0.13 



Some discussion and a couple of comments need to follow here. 

(i) In all cases, and particularly in the uj expansion approach, is decreasing 
very slowly. We have included the inaccessible value F = 25 in the table, just 
to illustrate this slow decrease. In the case of the NL in uj kernel, its value 
is practically constant in a certain Y region. For example, in the region of 
phenomenological interest 10^^ ^ 2; > 10~^, which corresponds to 5 < F < 9, 
we effectively have A^, = 0.30 -f- 0.29 and therefore the energy dependence 
of the saturation momentum is well described by the simple exponential law 
QliY) oc A^exp(Asl^). This is in agreement with saturation based models ||10 



that have successfully interpreted much of the HERA data. In these models 
such an ansatz for Qs was taken, with As as a free parameter, and after the 
fitting to the data was performed, the resulting value of A^ was very close to 
the one we quoted above. 

However, we should emphasize again that a fixed coupling description (which 
gives a constant value for A^, but much larger than the one found here) is not 
the right approach to the saturation problem, at least conceptually. The varia- 
tion of \s is rather slow partly because of the particular interplay between the 
leading behavior and the saturation corrections, as one can check in Eg. (|58|) , 
and partly because of NL corrections which slow down the evolution. 
Now we come to the issue of the actual value of Qs- As already mentioned 
at the end of Section 5, we do have control of the energy dependence of Qs, 
but there is an ambiguity of a free multiplicative constant as can be seen 
in Eg . (|53|) . Thus we cannot give any accurate prediction for the value of 
Qs- We just mention though, that for a "reasonable" choice of c in (|^), we 
obtain a phenomenologically accepted result. For example, with c = 3/4 (and 
A = 200 MeV), we find Q,(x ~ 3 x 10"^) ~ 1 GeV. 

In a very recent work a fit to the HERA data has been performed. The 
model relies on the combination of the Balitsky-Kovchegov and the DGLAP 
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equations. In such a case NL effects are neglected. As a part of their results, a 
best fit value of = 0.36±0.04 is given, where the form oc exp(AsF) has 
been assumed. It is interesting that the central value of this result coincides 
with our Leading in u one. 

(ii) It is clear that as Y increases the results of the u expansion converge to 
those of leading BFKL. This is expected as already mentioned at the end of 
Section 3. But this kind of convergence is slow, and it will not happen until Y 
reaches a very large value. Thus, in the region of phenomenological interest, 
the effects of the collinear improvement of the BFKL kernel are important. 
The leading BFKL prediction is reduced by an amount around 25 — 30% when 
NL in u corrections are included. 

It is also instructive to compare the leading and the NL in u results. In the last 
column of the table we give the ratio Sx/xo, which is the NL correction to the 
kernel divided by the leading one and evaluated at the point (7c, cUc). We notice 
that the correction has a "natural size" of the order of 0{as), when compared 
to the leading contribution, where the running is evaluated at the relevant 
value of Q^, say somewhere between Ql{Y) and Essentially this is a 

key element of the RG improvement of the BFKL equation. NL corrections 
are naturally small. We mention here that this observation is also true in the 



problem of the calculation of the hard Pomeron intercept cup ^ . 

(iii) In all cases, the second term in (|58|) which represents the effects of the 
saturation boundary, is smaller in magnitude than the leading term, but its 
contribution is extremely important. Even though pc is big, this term will not 
be negligible until Y, and therefore pc, reaches a huge value. 

It is quite interesting to investigate the parametric form of this saturation 
correction. In view of Eq.(p2D we naturally define = {I — 7o)Ac and lUs = 
(1 — 7o)As. Expressing Ac in terms of Pc in the 2nd term of (pBD, by keeping 
only the very leading behavior as determined from the dominant term in (|65|) , 
we obtain 



6^oXo(7o) 1 



46 



P< 



5/3 



1/3 



[as(Pc)] 



5/3 



(6^ 



— 5/3 

Contributions of the type as due to the running of the coupling have been 
observed and studied in the calculation of the hard Pomeron intercept 
331, BH RBI There is a striking similarity between the above equation and 
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Eq.(4.9) in |2J]. In our case the x function is evaluated at 70 and not at its 
minimum, and the additional 1/4 factor is understandable through Eg. (pT|) . 
We need to mention though, that in the saturation problem there is no limit 
in the applicability of ( p8D or ( |68D . On the contrary, since the line of evolution 
is one of increasing pc and decreasing ots-, the second term is always smaller 
than the leading one. 



7 The amplitude and geometrical scaling 

In this section we return to discuss the form of the amplitude. Since this 
is constant along the saturation line, we naturally expect a scaling behavior 
in a certain region around this line. Of course we restrict ourselves to the 
perturbative side of the boundary, since we do not know how to handle the 
details of the non-linear effects. The scaling form of the amplitude should be 
more and more exact as Y increases, but it is also a reliable expression even 
in the region of phenomenological interest. 

Returning to Eg. (^Sj) , with [3 and a as obtained in (pT]), and expressing the 
variables in terms of the more physical variables {Q,Qs), following the 
definitions we have given throughout this paper, it is straightforward to find 



T oc 



Q_ 



2 \ r 



Q\ l + c 



(69) 



Ql 1-70. 

The target dependent factors that are not written in this expression are of 
the form (const. //i^), with fi characterizing as always the target size. We 
notice also that the form of the amplitude will not change if we change Qg 
by a multiplicative constant, since such an operation can be absorbed in a 
redefinition of the constant c. Eg. (BUI) is valid in the region 



log§^< j^/^opy^-^, (70) 

a constraint arising from the expansion of the Airy function. Parametrically 
this constraint is satisfied since pc is growing with Y. In the region of interest, 
the number on the right hand side is not large enough to allow scaling in a 
wide region of values, as it has been observed in the HERA data [p7 |. 



Even at the very small value x ~ 10~^, the ineguality in ( |70D will totally break 
at ^ 50 GeV^ (independent of the unknown constant c). This "scaling" 
window depends relatively strongly on the initial condition. Furthermore, a 
more complete solution of ( ^9]) could also change the above estimate, but we do 
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not know how to control this issue in a more precise way. Other approximations 
that we have made, for example r] -C Pc, are also parametrically satisfied and 
are much weaker numerically. 

Eq.(|69|) is the same (but with a different Qs) as the one derived in Ref . |I8| , 
in a leading BFKL description of the saturation problem with either fixed or 
running coupling. Therefore it is interesting to try and interpret this result. 
The power behavior, with the anomalous dimension 1 — 70, is a consequence 
of the BFKL dynamics, either the "standard" or the collinearly improved. 
The logarithm in the square bracket is the effect of the absorptive boundary. 
Such a boundary imposes the condition that somewhere beyond the saturation 
line (where of course (^) is not valid any more since non-linear effects are 
important) the amplitude vanish. Clearly a logarithm satisfies this condition. 
We can give a more mathematical argument that ( |55D is the only scaling form 
allowed by the dynamics we have considered. 

To illustrate the issue in a simple way let us consider fixed coupling leading 
BFKL dynamics. In this case Eq.(^) simplifies to 



Following a similar procedure as in Section 3, we define now rj = p — Ps(V) 
and write the amplitude as T = exp[— (1 — 70)^7] "ipiv) to obtain 



with 7o a constant undetermined for the moment and where we dropped the 
d/dY term in order to look for an exact scaling solution. Clearly ps has to be 
constant. Expanding the Xo function around 70, (|72|) becomes 



Let us assume that that ip{ri) is a finite polynomial of rj and of order k, that is 




(71) 




(72) 




(73) 



i^iv) = 0-0 + a-iV + a2v'^ + ••• + o^kv'^- 



(74) 
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It is easy to see that the coefficients of the rf^ and r]^ ^ terms in ([73|) will vanish 
if and only if 

asXo(7o) - (1 - 7o)Ps = 0, (75) 

"sXo(7o) +Ps = 0, (76) 

respectively and 70 becomes our usual one. But if A; 7^ 0, 1 the coefficients of 
the remaining terms will not vanish and (^) is not satisfied. Thus only the 
case /c = 1, which is what we have found in (|69|), is allowed. The case = is 
allowed too, but it is not consistent with our saturation boundary conditions. 
Notice that ([75| ) and (^) are equivalent to the reduction of (^) and ( pTj ) in 
the fixed coupling BFKL case, with 7c — > 70 and — >• 0. 

The above can be generalized in the case of an RG improved kernel with 
running coupling in a straightforward way. Then the solution with /c = 1 is 
an approximate one when 77 -C and we are lead to Eqs. (^6[) and (|27|) , again 



with 7c — >■ 7o and /3 — >• 0. The one is satisfied in the asymptotic region and 
the other determines the leading behavior of p^. Notice that this procedure 
gives us a very fast way to obtain (|69D, but only the leading behavior of Qs, 
which is Qci is recovered. That is, we do not find in the corrections due 
to the boundary, which turned out to be important as we saw in the previous 
sections. 



8 Estimation of errors 

Here we check the approximations that we have made and which could be 
possible sources of errors. We find that all these errors are small and under 
control. 

(i) There is an ambiguity arising from the initial condition. Pc(0) should char- 
acterize the target size and is given by Eq . (p7|) . If we consider the scattering 
on a proton we may choose p to be around 1 GeV but we do not know the ex- 
act value. The same applies for A which can be taken to be around 200 MeV. 
However, changing Pc(0) by small amounts does not affect our results. We 
have already discussed that the corresponding corrections for are paramet- 
rically very small. From another point of view we ffist notice that Eq.(pO|) 
has "translational" invariance in Y . Thus, changing Pc(0), for example from 
2 log 5 to 2 log 3 which describes equally good a hadronic target, corresponds 
to a shift of the rapidity by a unit or so. Recalling our results from Section 6, 
where we found that is effectively constant in a wide regions of Y values, it 
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is clear the error is small. On the contrary, the value of Qs is sensitive to the 
initial conditions assumed. 

(ii) In Section 3 we have expanded our basic pc and 7c equations, ( ]26|) and 
(p7|), to first order in 57 = 7c — 7o- For the Leading in uj kernel, we easily find 
that 



^ 2Ac(l-7o)[l + Kl-7o)(l + Ac)] ^ Ac(1.85 + 0.591Ac) 
^ (l-7o)3x[,'(7o)(l + Ac)-2Ac 12 + lOAc ' ^ ^ 

If we expand to second order in ^7 we have to add on the right-hand side of 
Eq.(|0]) the term 

(57)V^'°(7o,Ac) 

26(1-70) • ^''^ 

This term gives a very small contribution in (^). Numerically it is 1 — 2% of 
the leading term when we plug in the expected value of Ac, the reason being 
the large numerical value of Xo(7o) which appears in the denominator in ([77|) . 
What is more interesting, is that this kind of correction can be neglected by 
making use of a more rigorous argument. From the equations above we see 
that the contribution is of order (!?(A^), which is of order O{uo'^) in view of 
(p^). Clearly this overlaps with NNL in uj corrections. 

(iii) In Section 5 we fixed the amplitude to be constant on the saturation 
line. This will be exactly true in the asymptotic region. In the intermediate 
region factors of the form exp[0(l/i(:^)] = exp[(yf(i:)] remain, because the func- 
tions hi and /12 are not equal to 1. These functions though, turn out to be 
close to 1, and therefore it is harmless to ignore these terms. Equivalently, 
we could multiply our solution (|43D with exp[— (yf(t)], which gives an equally 
good approximate solution since the remainder ( ^4|) changes just by a term 
hi{t)dg(t)/dt = 0{l/t'^). 

There are also certain assumptions and approximations that we made when 
deriving and ps as given in (^) and (^) respectively. We kept only the 
linear term in 5^ in the expansion of the Airy function in (|47|) , we replaced D by 
its asymptotic value Dq and we neglected ^^/t and ^/t^ terms in the exponent 
of the amplitude. If we relax these assumptions, we find that we should add 
in terms of the form 0{pc'^^^ , p~^ , ...). Notice that there are no 0{pc^^^) 
terms. All these corrections vanish when p^ becomes large. Thus in ( p2D all 
the terms that survive in the asymptotic region, including the constant one, 
have been kept. For completeness, we have checked that if we include all the 
available terms up to order C(p~^), then A^ is decreasing by ~ 0.01 in the 
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intermediate Y region. Therefore, we wouldn't like to sacrifice the relative 
simplicity of (|^ ) and (^Hj) for such an unimportant error. 

(iv) A crucial assumption that we made in order to simplify the BFKL equa- 
tion, was to drop the d/dY term that appeared inside the argument of the 
BFKL kernel in (^4D (of course we did that after switching to the scaling vari- 
able rj] the argument cannot be applied directly in ([17D ). One can understand 
that this term has a small effect, since we are looking for a scaling solution. 
The dominant d/dY term which appears on the left-hand side of (0) is al- 
ready one of the weak terms of our equation and therefore all other derivatives 
of this kind, which originate from NL corrections will be even weaker. It is 
interesting to understand this through a more rigorous argument. Let us con- 
sider the most important contribution of such a term, which is obtained by 
letting d/dr] —(1 — 70) and expanding to 1st order in the 2nd argument in 
(H). We get 

as a contribution to the left-hand side of our equation. For the case of the 
leading in u kernel, 0"'^ is 

We already see from ([79|) that the correction is suppressed by 1/pc, when 
compared to the dominant d/dY term. Then we just follow the procedure 
of Sections 3 and 4 and we find that ip{C,,t) is still given by (|43|) , where now 
the functions hi{t) and h2{t) are more complicated but their asymptotic limit 
remains the same. At the same time the pc equation changes because of the 
rid/ dY in ([79| ) . This will add on the right-hand side of Eq. (|30| ) a term 

'=0(7o,A,)-2feAe- ^ ^ 

However, as in case (ii) this is of order 0(c<j^) which overlaps with the NNL 
in u corrections. 

(v) We have considered three active flavors. A fourth quark, the charm, cannot 
be treated as massless as required by BFKL dynamics. Even if we set Nf = 4, 
then As will change by ^ 0.01. 



25 



9 Summary 



We have studied the energy dependence of the saturation momentum, using 
the colhnearly (RG) improved BFKL equation as our QCD dynamics. A sat- 
uration boundary of constant amphtude for dipole-hadron scattering has been 
used to describe the non-hnear effects in a simple but effective way. Follow- 
ing the evolution of the system, we were able to determine this boundary as 
far as the energy (or rapidity Y) is considered. The saturation momentum 
and its logarithmic derivative depend on F in a complicated way, par- 
ticularly at the Next to Leading Level. However, we have found that is 
decreasing very slowly, so that it is practically constant for a wide region of 
Y values and therefore is effectively described by the simple exponential 
law oc A^exp(Asy). In the region of current phenomenological interest 
we found at the NL level A^ = 0.30 -7- 0.29 with very reliable accuracy. The 
boundary (which imitates the non-linear effects), the RG improvement of the 
BFKL kernel at Leading and NL level, and the running of the coupling, were 
all important in obtaining Xs{Y). The exact value of Qg cannot be determined 
in our approach, because of a (not totally though) free multiplicative constant. 
The amplitude exhibits geometrical scaling in a certain region on the perturba- 
tive side of the boundary, and its form is independent of the particular BFKL 
dynamics considered. 
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Appendix 

Here we construct the NL in u kernel with an energy scale choice sq = Q 



following Refs.[22, 23, 24, Ml. We start with the NL BFKL kernel 19, gO 



Xiil) = - ^Ixlil) + X'oil)] - 7X0(7) 



vr^ cos(7r7) 



4(1 -27)sin^(7r7) _ 
67 

36 ~ 12 



4- 

3 + 



1 + 



Nf\ 2 + 37(1-7) 
Ar|J(3-27)(l + 27). 



+ ^C(3) 



5 



18 N. 



^)xo(7) + 



4 sin(7r7) 



$(7), (A.l) 
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where 



^(7) = . ■ . dx—— Li2(x), (A.2) 

6 sm(7r7) Jq 1 + x 

with the leading BFKL kernel Xo(7) and b as given in Eqs.(D and (|) in the 
text. The pole structure is different that the one given in because now we 
will include the Nf dependence. We need first to give the splitting functions 
and their Mellin transforms, which are 

W = n^T- + ^ + 4l-^)+bSil- z), (A.3) 

p,, = ^y+i^~^n (A.4) 
1 

A,,{u)= / dzz^Pggiu) -- 
Jo ^ 

= fe-T^ + :T^-T7^ - [^(2 + ^) - ^(1)], (A.5) 
1 + 2 + 00 6 + 00 



A,g{uj)= / dZZ'^P.giuj) 



Nf ( 1 



Define the function A^iuj) by 



(A.6) 



C 

At[uj) = Aggiuj) + -^A,g{uj), (A.7) 

with Cp = {N^ — i)/'^Nc. We can write the pole structure of Xiil) ^ 



^'^'^ 2(1-7)3 72 (1-7)2 

-f(- + + finite, (A.^ 



.7 1-7. 

where F, which vanishes when Nf = causing the simple poles to go away, is 



Nf 
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So far this refers to a symmetric energy scale sq = Qfi and in order to transform 
Xi{l) to the scale Sq = we need to add the piece — xo(7)Xo(7o)/2 [|T^ pOj . 
This will remove the triple pole at 7 = and will double the one at 7 = 1. 

Now, our leading in ut kernel xoi'y,^^) as given in (^) contains NL correc- 
tions which must be subtracted from our NL kernel to avoid double counting. 
Iterating u = agXoil,^) once, we easily find that we must add to the NL 
kernel the term xo{l)/{^ ~ lY- This will cancel the remaining triple pole at 
7 = 1. 

There are still double and single poles0. Putting all the previous pieces 
together and subtracting these remaining poles we find 



Xi(7) =Xi{l) - ^Xo(7)Xo(7) + ^^'^^^ 



1 /I 1 \ AriO) At(0)-6 



7 \7 1 — 7/ 7^ (1~7)^ 



(A. 10) 



a function with no poles at all at 7 = and 7 = 1. Of course we need now to 
add the correct form of all these subtracted terms that appear in the second 
line of the previous equation. We obtain 



Xih,^) =Xi{V + -- F[- + 



7 \7 1 — 7 + cij 
At{uj) At{uj) - b 
7^ (1-7 + tu 



Finally the full kernel is given as in Eq.(|T6| 
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